Hybrid Quantum Processors: 
molecular ensembles as quantum memory for solid state circuits 
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We investigate a hybrid quantum circuit where ensembles of cold polar molecules serve as long- 
lived quantum memories and optical interfaces for solid state quantum processors. The quantum 
memory realized by collective spin states (ensemble qubit) is coupled to a high-Q stripline cavity 
via microwave Raman processes. We show that for convenient trap-surface distances of a few 
jum, strong coupling between the cavity and ensemble qubit can be achieved. We discuss basic 
quantum information protocols, including a swap from the cavity photon bus to the molecular 
quantum memory, and a deterministic two qubit gate. Finally, we investigate coherence properties 
of molecular ensemble quantum bits. 
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During the last few years we have witnessed remark- 
able progress towards realization of quantum informa- 
tion processing in various physical systems. Highlights 
include quantum optical systems of trapped atoms and 
ions 1] , and cavity QED 2] , as well as solid state systems 
including Cooper Pair Boxes (CPB) ^ ^ and quantum 
dots 0. In particular, the strong coupling regime of cir- 
cuit CQED 6] was realized using a CPB strongly coupled 
to a strip line cavity. In the light of these developments 
it is timely to investigate hybrid devices with the goal 
of combining the advantages of various implementations, 
i.e. to build interfaces between, for example, a quantum 
optics and solid state qubit with compatible experimen- 
tal setups 7] . Such interfaces are particularly important 
in applications where long-term quantum memories or 
optical interconnects are required 
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In this Letter we study such a scenario by coupling a 
strip line cavity to a cloud of cold polar molecules |inl |. 
The cavity may be part of a solid state quantum pro- 
cessor involving CPB as charge qubits [3| and microwave 
photons as a quantum data bus, while the molecular en- 
semble serves as a quantum memory with a long coher- 
ence time. The condition of strong coupling between the 
cavity and the molecular cloud is achieved via the (large) 
electric dipolc moments of polar molecules for rotational 
excitations in the electronic and vibrational groundstate, 
which are in the tens of GHz regime, and thus provide 
an ideal match for resonance frequencies of strip line cav- 
ities. By adopting a molecular ensemble instead of a 
single polar molecule |ll| . we benefit from the enhance- 
ment of the coherent coupling gyN with the number 
of molecules N and g the single molecule vacuum Rabi 
frequency. This strong coupling between the molecular 
ensembles and the circuit CQED system also opens the 
possibility of a solid state based readout of molecular 
qubits. In addition, qubits stored in the molecular en- 



semble can be converted to "flying" optical qubits, using 
techniques demonstrated for atomic ensembles 12]. This 
provides a natural interface between mesoscopic quantum 
circuits and optical quantum communication. 

Let us consider the setup of Fig. ^ where two molec- 
ular ensembles are coupled to a superconducting cav- 
ity. The cavity is assumed to be strongly coupled to 
a CBP representing a circuit CQED system, as realized 
in recent experiments at Yale l6|. As discussed in de- 
tail below, molecular spectroscopy allows us to identify 
longiived states, for example in the form of a spin qubit 
|0), |1) in the ground rotational manifold. Starting with 
a cloud of N molecules prepared in |0)m = IO1O2 . . . Oat) 
coupling to a microwave or cavity field leads to exci- 
tations in the form of symmetric Dicke states, |l)m = 



1/v A^X^i |0i--li--0w) =^f^\^)m etc. For weak excitation 
the operator m obeys approximate harmonic oscillator 
commutation relations [m, m^] « 1, and the ensemble ex- 
citations are conveniently described as a set of harmonic 
oscillator states |0)m, |l)m = m^\0)m etc. Our goal be- 
low is to use the lowest two of these states as ensemble 
qubits, which can be manipulated by coupling them to 
the superconducting cavity and a CPB. 

The dynamics of the coupled system (Fig. la) can 
be described in terms of a Hamiltonian Hgys = He + 
Hm + HcM, which is the sum of a Jaynes-Cummings 
type Hamiltonian for the circuit CQED system He, a 
Hamiltonian for the (spin) excitations of the molecular 
ensembles Hm , and the coupling of the molecules to the 
cavity Hcm- The CQED Hamiltonian has the form 

He = -S,it)\e){e\+g,{\e){g\c+\g){e\c^). (1) 

Here \g) and |e) denote the ground and the first excited 
eigenstate of the CPB at the charge degeneracy point 
representing a charge qubit with a (tunable) transition 



frequency ujcq{t) with detuning Sc{t) 
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the cavity frequency ujc- The operator c (c^) is the cavity 
annihilation (creation) operator for microwave photons 
of frequency u}c^ and Qc is the vacuum Rabi frequency. 
Eq. ^ is written in an interaction picture with the bare 
cavity Haniiltonian UcC^c transformed away. The Hamil- 
tonian describing the internal excitations of the molecu- 
lar ensembles i — 1,2 and the coupling of ensemble states 
to the stripline cavity takes the form 



Hm + HcM 
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m^m, + ^5«(i)m|c + H.c. 



Before we enter the details of the derivation of Hm 
and HcM we note the basic structure of the Haniil- 
tonian Hsys- The ensemble excitations and the cavity 
represent a system of coupled harmonic oscillators in- 
teracting with a two-level system (CPB) with control- 
lable coefficients. In general, this provides the basic in- 
gredients for (i) swap operations between the charge, 
cavity and ensemble qubits, (ii) single qubit rotations 
of the molecular qubit via the charge qubit, and (iii) 
2-qubit entanglement operations between two ensemble 
qubits via the cavity mode, where the charge qubit plays 
the role of a nonlinearity For example, the CPB can 
act as a "single photon source" , i.e. we generate a 
superposition state of the charge qubit, which by an 
appropriate control sequence can be swapped over to 
the cavity, and is finally stored in one of the molecu- 
lar ensembles, and vice versa: {a\g) -\- f3\e)) |0)c|0),„ -^ 
\g) (a|0), + !3\l)c) |0)™ ^ |<?)|0), (a|0)™ + /3|1)™). 

The above discussion has ignored various sources of de- 
coherence. In the Yale experiment [g , the circuit CQED 
system realizes the strong coupling regime with vacuum 
Rabi frequency gc ^ Stt x 50 MHz. The decoherence 
of the charge qubit is dominated by the dephasing rate 
Tj"""^ « 27r X 0.5 MHz, while the photon loss rate is 
k/2'k = 1 to 0.01 MHz, i.e. the charge qubit is the dom- 
inant source of decoherence. Below we will show that 
for a cloud of iV « 10^ to 10^ molecules of temperature 
1 mK, trapped at a distance 10 /im above the strip line 
cavity, one can reach the regime of strong cavity - en- 
semble coupling 5m/27r « 1 to 10 MHz, which should be 
compared with the expected coUisional dephasing rates 
of the molecular memory of a few hundred Hz. 

Fig. lb shows the rotational spectrum of CaF, which 
provides an example for spectra of alkaline-earth mono- 
halogenides with a ^Ei/2 ground state corresponding to a 
single electron outside a closed shell. The spectrum con- 
sists of rotational eigenstates, described by a rigid rotor 
Hamiltonian Hr = BN^_^ with B -- 27r x 10 GHz the ro- 
tational constant, and N the angular momentum of the 
nuclei. The unpaired spin is coupled to the molecule rota- 
tion according to Hsr = jsrSN with •jsr '--^ 27r x 40 MHz 
and S the electron spin {S — 1/2). Coupled eigenstates 
are denoted by | A^, S, J; Mj) with J = N + S. As seen 
from Fig. ^ there is a spin rotation splitting (p-doubling) 
for rotationally excited states. In addition, there can be 
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FIG. 1: a) Ensembles of polar molecules and a CPB are cou- 
pled via the quantized field of a strip line cavity (see text for 
more details), b) Rotational excitation spectrum of molecules 
with a ^Ei/2 ground state where we plot excited states accord- 
ing to the Hamiltonian H = Hh + Hsr- For non-zero nuclear 
spin / (here I — 1/2) the hyperfine interaction leads to an 
additional splitting (& ~ 27r x 100 MHz) into eigenstates of 
F = J + I (hyperfine splitting for the excite states is not 
shown). Two qubit states in the N = manifold, |0), |1) are 
coupled by a Raman process involving a single cavity photon 
and an external microwave field, i}{t). 



hyperfine interactions, which in particular lead to a split- 
ting of the ground state A'' = 0, as in the case of CaF with 
a nuclear spin / = 1/2 which are coupled with J = 1/2 to 
F = and 1 states. In the following, we denote by |0), |1) 
a pair of states in the rotational ground state manifold to 
provide our spin qubit. Compared to qubits stored in the 
rotational degrees of freedom this choice of states avoids 
unfavourable N — 1 ^ N = collisions while |0) and |1) 
can still be coupled efficiently by a Raman process. 

The cavity mode and microwave fields of appropri- 
ate frequency and polarization couple rotational ground 
states to excited states with electric dipole matrix ele- 
ments /i (~ 5 Debye) . Two microwave driving fields pro- 
vide an effective coupling Hamiltonian ifJoff (0|0)(l|+h.c. 
to rotate the single molecule spin qubit, where fieff — 
r2ir^2/2A, with rii,2 the Rabi frequencies and A the 
detuning from the excited state |r), (A > rii,2). By 
similar arguments the coupling to the cavity has the 
form goff(t)|l)(0|c -I- h.c. with gcs{t) = gi^{t)/2A, where 
g = iiEc is the vacuum Rabi frequency, and Ec ~ 
y^huJc/2TreocPL is the electric field per photon for a cavity 
length L and typical electrode distance d. Typical values 
are g/2TT ~ 5 — 10 kHz for /z « 5 D and d « 10 /im. The 
distance d is also an estimate of the trapping distance 
of the molecular cloud from the cavity, which is well in 
the limit where standard trapping techniques work reli- 
ably and surface effects are negligible. Rewriting Hcm 
in terms of the collective operator m we obtain an effec- 
tive cavity - ensemble coupling gm{t) = VNg(.g(t). Due 
to the large wavelength Ac ~ 1.5 cm a trap volume of 
V < dx d X Ac/10 contains N = 10^ . . . 10^ molecules 



for gas densities of n ~ 10^^ cm -^ resulting in a coupling 
strength of gm/ 27r '--^ 1 - 10 MHz. The parameters Sm (t) 
are Raman detunings, which can be controlled indepen- 
dently, e.g. by applying local magnetic and/or electric 
fields. Thus we obtain the Hamiltonian Hm + Hcm, 



which allows a SWAP of a cavity and an ensemble state, 
for example, by an adiabatic sweep of Sm across the res- 
onance. This corresponds to a read / write operation 
Pc ® |0)m(0| < — > |0)c(0| ® Pm with pc an arbitrary den- 
sity operator of the microwave field in the cavity, and pm 
the identical state stored in ensemble excitations. 

The CPB provides a nonlinear element in the Hamilto- 
nian i/gys • This allows first of all single qubit operations 
of ensemble qubits, e.g. by combining a swap opera- 
tion with single qubit rotations of the charge qubit, and 
second, deterministic entanglement operations of qubits 
stored in two molecular ensembles. An example of such 
a protocol, which uses the CPB as a nonlinear phase 
shifter, is given as follows. We assume that the system is 
initially prepared in the state \ip)t=o = W m\0) c\g) with 
the charge qubit far detuned from the cavity resonance, 
1 5c (0)1 2> gc, and the two ensemble qubits in an arbitrary 
state \ip)m spanned by the basis |eie2)m, <^i = 0, 1. In a 
first step, in analogy to the single qubit swap, the state 
|V')m is (partially) transferred to the cavity. Assuming 
symmetric conditions, gm — gm and 8m {t) — 5m {t) 
it is convenient to rewrite the ensemble state in terms 
of the (anti)symmetric operators nis/a = ("^i i ?n2)/v2 
acting on |00)m. An adiabatic sweep of the Raman de- 
tunings then realizes the swap operation m||00)m|0)c -^ 
\m)m\l)c and (mt)2|00)„|0)c ^ V2|00)™|2)c while the 
states |00),„|0)e, mt|00)™|0)e and {mif\m)m\Q)c re- 
main unaffected. In a second step, the charge qubit 
is adiabatically tuned close to resonance for a time T, 
1^0(2^/2)1 < gc- During this pulse the non-vacuum states 
acquire a nonlinear dynamical phase, \n)c — > e"''"|n)c 
with 4>n = ~j^dt'{5c{t') + ^5l{t')^nAgl)/2. The 
pulse form 5c{t) and the length T are chosen such that 
</>! ~ 7r/2 and 02 = 27rn (see e.g. Fig. EJ. The 
second condition ensures that after writing the cav- 
ity state back into the ensembles states, i.e. reversing 
step one, the ensemble states |20)m and |02)m remain 
unpopulated. The total gate sequence corresponds to 
a 'VS WAP" -like gate for two ensembles qubits, with 
|00)^ ^ |00)™, |10),„ ^ eW4(|io)„^ + i|01>„0/V2, 
|01)„ ^ e'^'\i\lQ)m + |01),„)/V2 and |11),„ ^ |11)™. 
A numerical simulation of this gate sequence based on 
a master equation treatment of the dissipative terms |6| 
shows that the gate fidelity is only limited by {gcT2)~^ 
i.e. the decoherence of the CPB during the time it is 
tuned close to resonance (see Fig. |21for more details). 

We now turn to an analysis of decoherence in the 
molecular ensemble. In particular coUisional dephasing of 
the ensemble qubit, and a spatial variation of the cavity- 
molecule coupling geS_{x) in combination with the ther- 
mal motion of the molecules in the trap contribute to a 
finite decoherence time of the molecular quantum mem- 
ory, and result in imperfections during gate operations. 

An operational definition of the decoherence time of 
the ensemble qubit can be given in terms of an (ideal- 
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FIG. 2: a) Adiabatic energies levels for the states \n)c\g) as 
a function of the charge qubit detuning 5c- b) Evolution of 
the dynamical phases (j)\ (solid line) and cj>2 (dashed line) for 
a pulse 5c{t) = -<5o(2t/r - 1)^ - 5i , Sq/qc = 30, Si/g^ = 0.44 
and T = 44.79/(?c. For the same pulse 5c{t) the resulting 
fidelity of the total gate sequence, Ta (averaged over all initial 
states \tp)m), is plotted in c) for diflterent values of the charge 
qubit dephasing rate, Tj"^, and the cavity loss rate k. 



ized) experiment. A cavity qubit Pc{t = 0) = |'0)c('!/'| 
with |-0)c — a|0)c 4- /?|l)c is written at time t = to the 
molecular memory with all molecules initialized in the 
state |0) in a (perfect) swap operation, kept in storage 
for a time interval r, undergoing dephasing collisions. At 
t = T, the qubit is transferred back to the cavity mode, 
resulting in a reduced density matrix Pc{t) of the cavity 
with a fidelity _FV — min^^ c{4'\Pc(t)\'iI))c, with the deco- 
herence time of the ensemble memory identified as the 
decay time of the fidelity. The analysis of this process re- 
sembles the discussions of clock shifts, and in particular 
studies of collisional dephasing of spins in thermal and 
quantum degenerate atomic clouds in a Ramsey intcrfcr- 
ometry setup. A formal theoretical description of these 
phenomena is provided by quantum kinetic theory Il3l |. 
We consider a thermal cloud of molecules in the lowest 
rotational state with qubits stored in spin or hyperfine 
states (Fig. ^ . Molecules in the rotational ground state 
are trapped magnetically [lj| or by a (spin independent) 
electric rf-trap 15] . In a magnetic trap two molecules 
interact asymptotically according to a I^(r) ~ —Cq/t^ 
potential with Cg = (/i^/47reo)^/6i3. The effective range 
of this potential is given by i?* = y^mCe/h?, which pro- 
vides an estimate for the s-wave scattering length a. For 
example, for CaCl (which has two magnetically trapped 
hyperfine states) we obtain R^ w 780 as which is a few 
times the typical scattering length encountered for al- 
kali atoms. S-wave scattering dominates for tempera- 
tures T < T, « 1/iK where the thermal energy is below 
the centrifugal barrier for higher angular momenta, lead- 
ing to an estimate for the collision rate 7coi = Sna'^nv « 
2tt X 150 Hz for n = lO^^cm"'^ and v the relative thermal 
velocity. For temperatures T ^ T^ higher partial waves 
contribute, and an estimate of the cross section based on 
the unitarity limit gives 7coi ^ 27rx 700 Hz for T = 1 mK. 
In electric traps the induced dipole moments, /iind, lead 
to a V{r) ^ A'fnd/''^ dependence of the asymptotic inter- 
action. Although this long-range behavior significantly 
changes the low temperature scattering (T < 1/iK) the 



estimate based on the unitarity limit at T sa 1 mK still 
provides a valid bound for scattering rates of weakly po- 
larized molecules (/lind < ID). 

We have calculated the decoherence time of an en- 
semble qubit corresponding to coUisional dephasing using 
quantum kinetic theory 16]. Dephasing of the qubit co- 
herence {pc)io{t) = exp(— 7iot/2)(pc)io(0) is associated 
with spin dependent collisions between the states |0) and 
|1), which gives a contribution 
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X (|/o^o(K) - /o^(K)p + |/™(K)p + |/^?(K)n . 

It depends on the difference between /qq and f^i, the 
elastic scattering amplitudes for the internal states |00) 
and (1 10) + |01))/-\/2 averaged over the thermal distri- 
butions P{k) in the scattering process between momenta 
K = (^3,^4 <— ^1,^2), and i5- functions accounting for 
energy and momentum conservation in the collision. In 
addition, there may be contributions from inelastic col- 
lisions, /oo and /qq, which scatter molecules outside 
the |0), |1) subspace. While accurate scattering ampli- 
tudes for molecular collisions may not be available at 
present, we can estimate these contributions in certain 
limits. For s-wave scattering the above expression sim- 
plifies to 710 = 87r(aoo — aoi)^"'^ with aoo and oqi scat- 
tering lengths. If we assume that the scattering length 
is dominated by a spin exchange potential, the scatter- 
ing is characterized by a singlet (as) and triplet scat- 
tering length (ay). In the simple case of a pure spin 
qubit {|0),|1)} = {|5 = l/2,ms = ±l/2)} we find aoo = 
Qoi = clt, and the dephasing rate is determined by non- 
vanishing contributions arising from magnetic dipole and 
spin rotation coupling, which are expected to be much 
smaller. In a similar way, in the presence of hyperfine in- 
teractions we can form a qubit |0) = |F = /-|-l/2, Mp = F) 
and |1) = |F' = /-1/2,Mf = F'), where again |00) and 
|01) -|- 1 10) contain no spin singlet contribution and the 
leading dephasing term vanishes. In the worst case the 
decoherence rate is bounded by the single molecule colli- 
sion rate 710 ~ 7coi which has been estimated above. 

Spatial variations of the effective single molecule-cavity 
coupling, 3eff{a;), result in a dephasing of the qubit 
during a single swap gate and an incomplete recovery 
of the state after a redistribution of the molecules be- 
tween two successive write/read operations. The in- 
homogeneity in the coupling arises from the variation 
of the cavity mode function on a scale of the elec- 
trode distance, g{x) « g(\ — ax/d)^ with a a numeri- 
cal constant, and a position dependence of the detun- 
ing A(a;) « A — mSuj'^x'^ /{2H). Here (5w^ = ujf — lu^ 
accounts for a difference in the trapping potentials for 
the qubit states (ujt) and the excited state \r) (ujr)- For 
an optimal detuning A* ~ ^ig"^ N [kiJT Suj'^Y / Hoo^ti^ the 



inhomogeneous coupling results in a total gate error of 
e « a^ikbT/mco^d^) + (kbTStu^K/hg^Nuj^)^/^ 1^. For 
gVN - 27r X 10 MHz, k - 27r x 10 kHz and at T = 1 
mK gate fidelities oi T > 0.99 require trap frequencies 
of Wi '^ 27r x 50 kHz and a similar trapping potential for 
the state \r) {Scu'^ ~ O.lw^). Lower temperatures and an 
optimized cavity / trap design, e.g. a, Suj'^ -^ 0, lead to 
a further significant reduction of gate errors. 

In conclusion, thermal ensembles of cold polar 
molecules represent a good quantum memory that can 
be strongly coupled to strip line cavities with long life 
time, limited essentially only by coUisional dephasing. 
We note that these dephasing channels can be virtually 
eliminated, if the ensemble is prepared in a high-density 
crystalline phase of dipolar gases with dipole moments 
induced and aligned by a DC electric field under 2D trap- 
ping conditions jig . The present work opens an exciting 
avenue towards long lived molecular quantum memory 
for solid state quantum processors. 
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